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We review what is presently known about higher loop corrections to the Euler-Heisenberg 
Lagrangian and its Scalar QED analogue. The use of those corrections as a tool for the 
study of the properties of the QED perturbation series is outlined. As a further step in a 
long-term effort to prove or disprove the convergence of the N photon amplitudes in the 
quenched approximation, we present a parameter integral representation of the three- 
loop Euler-Heisenberg Lagrangian in 1+1 dimensional QED, obtained in the worldline 
formalism. 
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1. Introduction 

Heisenberg and Euler's 1936 calculation of the one-loop effective Lagrangian in- 
duced for a constant Maxwell field by a spinor loop was not only a milestone in 
the development of QED, but remains until today the prototypical example for the 
concept of integrating out degrees of freedom in field theory. In practical terms, 
it encodes in a very concise form information on a host of photonic low-energy 
processes (see f21 for a review) . Higher loop corrections to the Euler-Heisenberg La- 
grangian and its Scalar QED analogue, due to Weisskopf ^ (both called "EHL" in 
the following for simplicity) were studied only much later, starting with Ritus' 1975 
calculation'^ of the two-loop EHL. The purpose of this talk is to give a summary 
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on what is known about the EHL at the multiloop level, and to argue that those 
multiloop corrections, although not likely to be of phenomenological interest in the 
near future, contain important structural information on QED. 

2. QED in a constant external field in the worldline formalism 

We start with a short introduction to the worldline representation of the QED 
S-matrix, going back to Feynman which in recent years has emerged as an 
extremely efficient tool for the computation of processes involving constant external 
fields in QED. For the simplest case of the one- loop effective action in Scalar QED, 
it reads'^ 



T[A] 



dT 

( 

T 



Vx exp 



dr ( + ieA^i^ 



(2.1) 



Here m and T are the mass and proper-time of the scalar particle in the loop, 
and the path integral runs over the space of closed trajectories with period T, 
x^{T) — x'^{0), in (euclidean) spacetime. The spinor QED equivalent of (2.1) is 



obtained ^ by the addition of a global factor of — and the insertion of a spin 
factor S[x, A] under the path integral. The modern way of writing this spin factor 
is in terms of an additional Grassmann path integral'^, 



S[x,A] 



ViPIt) exp 



(2.2) 



Here the path integration is over the space of anticommuting functions antiperiodic 
in proper-time, V'^(ri)V'''(r2) = -V{t2)'iP'"{ti), V{T) = -^'''(0). 

Presently, three quite different methods are available for computing such path 
integrals (see |H! for a review) : (i) the "string-inspired approach" , based on a per- 
turbative expansion and gaussian path integration EEoEiEEIl (ii) the "worldline 
instanton approach" , using a stationary path approximation EEII and (iii) the di- 
rect numerical calculation using Monte Carlo techniques 

All three methods have been applied to the calculation of Euler-Heisenberg 
Lagrangians. We will explain here only the "string-inspired" method; see IT^ IT?l [T71 
for the worldline instanton and [121 IIHI for the worldline Monte Carlo approach. If 
we expand the interaction exponential, 



exp 



dr ieA^x^ 



E 



(-ze) 



N N 



n 



dTi 



x^'{n)A^{x{T,)) 



(2.3) 



the individual terms correspond to Feynman diagrams describing a fixed number of 
interactions of the scalar loop with the external field, see fig. [TJ 

The corresponding N - photon scattering amplitude is then obtained by special- 
izing to a background consisting of a sum of plane waves with definite polarizations. 
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Fig. 1. Feynman diagrams equivalent to the one- loop effective action. 

— '^iLi^ifj.^^''''^ 1 ^i^d picking out the term containing every once. This 
yields the following representation of the N - photon amplitude, 



'dT 

I 

T 



(2.4) 



Here V^^-^ [k, e] = dr e • i(T) e^^^i"^) denotes the same photon vertex operator used 
in string perturbation theory. The integral over the zero mode XQ = }p J^^ dr (t) 
factors out and produces the usual energy-momentum conservation factor. The re- 
duced path integral / Vy^r) over 2/(r) = x{t) — xo is gaussian, and can be evaluated 
using the "bosonic" worldline Green's function Gb, 



{y^{ri)y''{T2)) = -5^^GB(ri,r2) = -g^'' | n - | 



T 



(2.5) 



Using a formal exponentiation of the factors £i ■ Xi 's and "completing the square" 
yields the following closed-form expression for the one-loop N - photon amplitude: 



r°° dT ^ 
T[{h,e.}] = {~ief{2^f5{Y,h) (47rT)-^e-™'^n / 

Jo ^ ,,^1 Jo 

r ^ 1 . 1 .. 1 

X exp<^ ^ [i^'^Bijki ■ kj + iGsijh ■ ej + -GBijSi ■ £j\ > 



multi— linear 



(2.6) 

Here it is understood that only the terms linear in all the £i, . . . ,£Ar have to be 
taken. Dots generally denote a derivative acting on the first variable, and we ab- 
breviate Gsij = GB{Ti,Tj) etc. The factor {AttT) ^ represents the free Gaussian 



path integral determinant. The expression (2.6) is identical with the "Bern-Kosower 



Master Formula" for the N photon case EH 



In the spinor QED case, the correlator for the evaluation of the additional 
Grassmann path integral is {^{ti)iP{t2)) — ^g'^"'Gp{Ti,T2), with Gf{ti,T2) = 
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sign(Ti — T2). Its explicit evaluation can, however, be circumvented, using the fol- 
lowing "replacement rule" Writing out the exponential in the master formula 



eq.(2.6) for a fixed number N of photons, one obtains an integrand 



exp 



{-i)^ PN{GBij,GBij) exp 



I multi— linear 



1 ^ 

7, X! '^Bijh ■ k.j 



(2.7) 



with a certain polynomial Pm depending on the various Gsij^s, Gsy's, as well as on 
the kinematic invariants. By suitable partial integrations all second derivatives Gsij 
appearing in Pj\j can be removed, so that Pn gets replaced by another polynomial 
Qn depending solely on the Gsy's, 



PwiGsij, GBij)e^ 



QNiGsij) e 



5 Yl Gsijki-kj 



(2., 



Then the integrand for the spinor loop case can be obtained by simultaneously 
replacing every closed cycle G Biii2G Bi2ii ' ■ ■ Gsi^ii appearing in Qat by 



G Bi\i2 G Bi2i3 



G 



Bik'ii 



G 



Fikii 



(2.9) 



An additional background field A^^ (x) with constant field strength tensor can 
be simply taken into account by appropriate changes of the worldline propagators 
and the path integral determinant. Those areESGlESEIl (deleting the "bar") 



G 



B12 



Gbi 



T 



'iZGBi 



2{zf Vsin(Z)~ 



iZG 



B12 



-iZGE 



Gf12 Gf12 — Gf12- 

{inPy^ -> (47rT)"*det" 
{'inPy^ -> (47rT)"*det" 



cos(Z) 
sin(Z) 
Z 

tan(Z) 



(Scalar QED) 
(Spinor QED) 



(2.10) 

These expressions should be understood as power series in the matrix Z^'^ = eTF^'^. 



Thus one obtains the following generalization of (2.6), representing the scalar 
QED N - photon scattering amplitude in a constant field [111121 



'dT 
~T 



^det"5 



sin(Z) 



N T r ^ rl 

/ dr^exp^ ^ ^-ki ■ QBij ■ 



kj — iSi ■ QBij ■ kj + -Ei ■ QBij ■ £j 



I multi— linear 



(2.11) 
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The cycle replacement rule (2.9) can also be generalized to the constant field case. 



The master formula (2.111 is valid off-shell, and can therefore be used to con- 



struct the quenched (one-electron-loop) Euler-Heisenberg Lagrangians in scalar or 
spinor QED at the n - loop order by starting from the one-loop 2(n — 1) photon am- 
plitude in a constant field, and sewing off pairs of legs. (Alternatively, one can use 
the worldline formalism also to calculate the Euler-Heisenberg Lagrangians directly 
at the muhiloop level IMHEl) 



3. The one-loop Euler-Heisenberg Lagrangians 



The one-loop EHL's involve only the determinant factors in (2.10). After renormal- 
ization, one has 



£^^^ (F) = - 

spinV / 



1 

167r2 
1 



dT 



spin 



T3 
' dT 



6 



(a^ - b'^)T^ - 1 



3 ^ 



- b^T^ - 1 



{eaT){ebT) 
sinh(ear) sin(e6T) 

(eaT)(e6T) 
tanh(eaT) tan(e6r) 

(3.1) 

Here a, b are the two invariants of the Maxwell field, related to E, 'B hy — b^ — 
B"^ — E^, a6 = E • B. The subtraction terms in the square brackets implement the 
renormalization of charge and vacuum energy. The subscripts distinguish between 
Scalar and Spinor QED, the superscripts refer to the loop order. 

The EHL's contain the information on the N photon amplitudes in the low 
energy limit where all photon energies are small compared to the electron mass, 
uJi <^ m. At the one- loop four photon level, this construction of the low-energy 
amplitude from the effective Lagrangian is a textbook exercise (see, e.g.,[53|. Even 
the one-loop (on-shell) N - photon amplitudes in this limit can still be written quite 
concisely using spinor helicity techniques 



Except for the purely magnetic field case, the parameter integrals in ( 3.1 ) contain 



poles, leading to an imaginary part of the EHL's. A simple application of the residue 
theorem gives Schwinger's famous formulas 



fc=i 



167r3 



fc=i 



cxp 



(-1) 



k + l 



fc2 



tt/c 

T 

exp 



ttA: 

T 



(3.2) 

(/? = eE/m?). The kth term in these sums is interpreted as representing the insta- 
bility of the vacuum with respect to the coherent production of k electron-positron 
(reps, scalar-antiscalar) pairs by the field (vacuum tunneling). In the following we 
will concentrate on the weak field limit /3 <C 1 where only the k = I term is relevant. 
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4. The two-loop Euler-Heisenberg Lagrangian 

The two-loop EHL, involving one internal photon exchange in the loop, was first 
calculated by V.I. Ritus, both for Spinor and Scalar QED This resulted in a 
type of rather intractable two-parameter integrals, on which also later recalculations 
were not able to substantially improve P ''|20|22 [ jjQ-^^gyer, the first few coefficients of 
the weak-field expansions of the two- loop EHL's have been computed pop^P^ ^ and 
there are simple closed-form expressions for the case of a (euclidean) self-dual field 
Those allow one to extend the one-loop calculation of the on-shell low energy N 
photon amplitudes, mentioned above, to the two-loop level for the case where all 
the photon helicities are the samel^. 



As to the imaginary parts, the Schwinger formulas (3.2 1 generalize to the two- 
loop level as follows fSiJ 



Im£S„(£;) = ^/3^ J^aTTiff exp 



k=l 

4 oo 



Trfc 



k=l 



Trfc 



(4.1) 



(a = 1;^) where 



seal, spin 



Ck 

'Vl3 



+ l + 0(v//3) 



k-l 



Ck — — ^ — , , fc > 2 



(4.2) 



Thus at two-loop the fcth Schwinger-exponential appears with a prefactor which is 
still a function of the field strength, and of which presently only the lowest order 
terms in the weak-field expansion are known. Still, things become very simple at 
leading order in this expansion: Adding the one-loop and two-loop EHL's, one finds, 
e.g. for the spinor QED case 



87r3 



(4.3) 



and this result is spin-independent (but for the normalization). In it was fur- 
ther noted that, if one assumes that in this weak- field approximation higher order 
corrections just lead to an exponentiation. 



^Im£(')(i?) Im£(i)(£;)e"^ 



(4.4) 
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then the result can, in the tunnehng picture, be related to the fact that a created 
pair gets born at a finite distance, and thus with a negative Coulomb binding energy. 

5. An all-loop conjecture from worldline instantons 

Already in earlier work (but unknown to the authors of [31] at the time) Affleck et al. 
had invented the concept of worldline instantons, and their principal application 



was precisely to demonstrate the exponentiation property (4.4), although not for 
spinor QED but for the scalar QED case. Even though neither derivation of (4.4) 
(to be called "AAM conjecture" in the following) can be considered rigorous, the 
fact that it was obtained by two very different lines of reasoning makes us confident 



about its correctness. Assuming this to be the case, (4.4 1 is a highly remarkable 
result, since, despite of the simplicity of the derivation by [14] in terms of a single 
semi-classical instanton trajectory, it is a true all-loop result, receiving contributions 
from an infinite set of Feynman diagrams (although it is important for the following 
that only quenched (one fermion-loop) diagrams contribute in this limit). Moreover, 



[TJ] argue that (4.4) is written in terms of the physically renormalized mass, that 
is, the counter diagrams for mass renormalization, which normally are necessary 
starting from the two-loop level, have been taken care of implicitly. This again 
cannot be considered as rigorously proven, however what can be shown easily is that, 



if (4.4) holds, then the mass appearing in it must be the physically renormalized 
one. Namely, using a Borel dispersion relation one can show^^S the following general 
relation between the prefactor of the first Schwinger exponential and the leading 
asymptotic growth of the weak field expansion coefficients at fixed loop order I: 
Defining these coefficients by 

oo 171 9 

£(')(£;) = ^c«(.)(^)" (5.1) 



one has 



c(')(n)"^-c«r[2n-2] 

(5.2) 

with constants c^- This implies, in particular, that the leading factorial growth 
order of the expansion coefficients must be the same at each loop order, and in [28] 
it has been shown that already at the two-loop level this holds true if and only if 
the renormalized mass is the on-shell one. 

6. Relation to the multiloop N photon amplitudes at large N 

Now, the AAM conjecture is remarkable not only for its simplicity, but also for 
the fact that, despite of arising from a true all-order loop summation, the result is 



December 6, 2011 1:22 WSPC/INSTRUCTION FILE qfextllarxive 



8 Authors' Names 

analytic in a. This appears to run contrary to many arguments which have been 
given, starting with Dyson's classic 1952 paper to show that S-Matrix elements 
in QED can never be analytic in a. Now the Schwinger pair creation rate is not 
itself an S-matrix element, but it can, at any loop order, be related to the N - 



photon amphtudes at large N using the above relations (5.1 1,(5.2) and the already 
mentioned standard procedure for converting the weak field expansion coefhcients 
into low-energy photon amplitudes (this involves also an extension of the AAM 
conjecture from the electric field case to the general constant field or at least self-dual 
field case 33^ ^^^i as other modest assumptions). This led one of the authors and 
G.V. Dunne to conjecture in 20041^ that the perturbation series for the N photon 
amplitudes, albeit divergent for the full amplitude, is indeed convergent at the level 
of the quenched approximation. This conjecture is not in contradiction with existing 
general theorems on the QED perturbation series, and extends a corresponding 
conjecture made by Cvitanovic in 1977 for the g — 2 factor. If true, it would 
indicate extensive cancellations between Feynman diagrams, presumably due to 
gauge invariance. 

Referring for the details to |34|30|33|i7]^ ^et us just state here that this line of 
reasoning leads to a number of nontrivial predictions starting from the three-loop 
level. Namely, we expect to find the expansion coefficients of the three-loop EHL 
(for both Scalar and Spinor QED) to have the following three properties: 



A3) 



(1) lim ^ (") - Iq.2 

(2) Only the quenched part of the EHL should contribute to this limit. 



(3) The convergence of to (from the first eq. in (5.2)) should not be 



slower than the one for the corresponding two-loop to one- loop ratio. 

Unfortunately, a calculation of the three-loop EHL has so far been proven tech- 
nically out of reach. However, motivated by work of Dunne and Krasnansky IHIIZ] 
on EHL's in various space-time dimensions it was shown in [T7] that the whole above 
machinery can, mutatis mutandis, be transferred to the computationally simpler 
context of QED in 1+1 dimensions. In particular, all of the above three-loop pre- 
dictions can be generalized to this case, changing only the definition of the ratios. 



c^^Hn) 4^2(n + /-l) 

2 

and changing a to 2^^. Preliminary results on a calculation of the three- loop EHL 
in 2D Spinor QED using the Feynman diagram approach were presented at the 
QFEXT09 conference however this approach ultimately failed, since it led to 
parameter integrals with spurious IR divergences. Here, we will present the results 
of a new run on the calculation of the same three-loop EHL using the worldline 
formalism!^, which has allowed us to obtain this Lagrangian in terms of manifestly 
(IR and UV) finite parameter integrals. 
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7. The three-loop Euler-Heisenberg Lagrangian for 2D QED 

At the three-loop level, there are three Feynman diagrams contributing to the EHL 
in spinor QED, shown in fig. 2: 





A 



B 



C 



Fig. 2. Feynman diagram representation of the three-loop EHL. 



Here the solid line stands for the fermion propagator in the constant field with 
field strength tensor F = ( q ) . The calculation of the non-quenched diagram C 
is straightforward, and leads to a fourfold proper-time integral 



167r3/ 



dzdz'dzdz" 



sinli z sinh z' sinh z sinh z" 
[sa\\\{z + z') sinh(£ -|- z")f 

2k(z+z' +z+z") 



sinh z sinh z' sinh(i + z") + sinh z sinh z" sinh(2: -|- z') 



(7.1) 

where k = m? /2ef. From this we have obtained, by numerical integration using 

MATHEMATICA, the first 12 weak field expansion coefficients, which was sufiicient 
to verify that they are indeed asymptotically suppressed, even exponentially, with 
respect to the asymptotic prediction of the (2D analogue of) the AAM formula. 
Thus point (2) of our three-loop predictions above has been settled. 

To the contrary, the parameter integral representations which wc have obtained 
for the quenched diagrams A and B are lengthy, and cannot be given here in full. 
Their structure is 



'dT 



z 



n 



t&nhZ -.-Jo 



-T 



dn 



X (2/1234 + /l324 + 4/123 + 2/12 + 4/13 + /i2,34 + 2/13,24) 



(7.2) 
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where Z = e/T and, for example, 



^ijkl — 



tr{{ijkl}i 



(7.3) 



with (compare (2.9)) 



{iii2 ■ . . 



^ Biii2 GBi2i3 



■ GF^„^, (7.4) 



Here we have further introduced Gsij ■— Gsij — Gbu + GFii, and A is a determinant 
also involving the worldline Green's functions. Note that (7.2) represents the sum 
of both diagrams A and B. 



8. Conclusions 

We have summarized here what is presently known about multiloop corrections to 
the EHL, concentrating on the potential of such corrections to yield information 
on the high-order behavior of the perturbation series in QED. As part of a long- 
term effort to prove or disprove "quenched convergence" for the case of the photon 
S-matrix, we have presented a parameter integral representation of the three-loop 
EHL in 2D QED. Although the extraction of the weak- field expansion coefficients 
from this representation and verification of the remaining predictions implied by 
this conjecture at the three-loop level (points (1) and (3) above) will still require 
very substantial work, we are confident that we will have definite results to show 
by the time of QFEXT13! 
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